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STRESS  ANALYSIS  ABD  STATIC  TEST  OF  A  THREE  CELL  BEAM 

IHTROPqCTION 

An  aircraft  wing  conalsts  of  structural  members  capable  of  realstlng 
bending,  torsional  and  shear  loads.  Insofar  as  the  stress  analysis  of  the 
ving  is  concerned,  the  structural  members  resisting  each  of  these  loads 
can  be  considered  separately.  In  the  consideration  of  bending  loads,  ving 
structures  may  be  typed  according  to  the  bendlng-load  resistant  material. 
For  example,  all  beading  material  may  be  concentrated  in  the  spar  caps  or 
dlBtril>uted  around  the  periphery  of  the  profile.  The  analysis  of  the 
first  type  treats  the  ving  covering  as  ineffective  in  bending.  This  is  a 
reasonaljle  conclusion,  inasouch  as  the  sheet-metal  covering  vill  buc]s:le  at 
a  very  lev  load  and  the  load  carrying  ability  insofar  as  beading  la  con- 
cerned, is  therefore  negligible   As  for  the  second  type,  the  distrll3Uted 
bending  material  consists  of  spanvise  stiffening  members.  These  members 
are  in  the  form  of  stringers,  corrugations,  and  inner  skins  fastened  to  the 
covering  metal.  Hovever,  on  the  compression  side  of  the  ving  the  metal 
sheet  covering  vlU  buckle  at  relatively  lov  loads  vhlch  necessitates  the 
use  of  effective  widths.  Hence,  the  stress  analysis  of  such  a  structure 
considers  the  skin  as  partially  effective  in  bending. 

At  present  there  is  the  necessity  of  oox^toructing  aircraft  vings  vith 
comparatively  thick  metal  sheet  covering.  This  increase  of  covering  thick- 
ness pend.ts  the  ving  to  carry  relatively  high  airloads  before  buckling 
of  the  skin  on  the  compression  side,  thus  the  skin  is  fully  effective  in 
bending.  The  stress  analysis  of  this  type  of  ving  is  the  subject  of  this 
report. 
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THREE  CELL  BEAI.4  WITH  SKIN  FULLY  EFFECIVE  IN  BENDING 


PART  I 

ANALYSIS 

THREE  CELL  BEAM  WITH  SKIN  FULLY  EFFECTIVE  IN  BENDING 


Tho  object  of  thle  analysis  Is  to  doTolop  tho  shear  flov  equations 
for  a  three  cell  beam  irlth  skin  fully  effectlTe  In  bending.  The  re- 
sulting equations  vlll  apply  for  a  beam  loaded  with  shear,  bending  and 
torsional  loads. 

Consider  a  three  cell  beaa  with  a  general  profile  (e.g.,  any  alirfoll 
shape)  as  shown  in  Figure  1. 


>'X 


Figure  1 

Beferlng  to  Figure  1. 

0  =  centrold  of  cross  section, 
z     =  distance  from  centrold  to  shear  center. 

B.C. 

V  =  shear  force  In  y  direction  (through  sheea"  center). 
Initial  Assumptions: 

1.  Principal  axes  are  y  and  x  axes,  where  x  axis  is  chord  line. 

2.  Wing  planform  contains  no  taper. 

5.  Load  applied  through  shear  center  in  y  direction. 
The  shear  flow  designatioiw  are  as  shown  in  Figure  2. 


3.3'  _^« 


Figure  2 


Referring  to  Figure  2. 


<1  =  Shear  flow  at  any  point  noae  skin. 

q^  =  Shear  flov  at  any  point  on  upper  skin  of  center  cell- 

a  =  Shear  flow  at  any  point  on  rear  cell  skin. 

a     -   Shear  flow  at  any  point  on  "bottom  skia  of  center  cell, 
a   =  Shear  flow  at  any  point  on  front  web. 
q^   =  Shear  flow  at  any  point  on  second  web. 

q  =  Value  of  a     at  point  a. 
q  ,  =  Value  of  q,^  at  x>olnt  a*. 

q  =  Value  of  q^  at  point  d. 

Cenaider  now  the  span  wise  forces  acting  on  section  p  a  of  Figure  2, 
as  shown  in  Figure  3.  Fron  equlllbriusi  of  spanwise  forces: 


dP 

dZ 


P  +  ~  -  P  +  q^dZ 


Figure  5 


q  dZ  =  0 
a 


\-\ 


dP 
dZ 


ap 


where  subscript  ap  indicates  re- 
gion from  point  a  to  point  b. 


-M  J 

f  =  -^ 
b    I 

z 


Thua:  P  = 


vbere  f,  ^  Stress  due  to  beollng. 

0 

J  -  Distanc*  from  x  axis. 
I  =  MoBont  of  Inartla  about  x  axis. 

X 

M     =  Bending  aoaent  about  x  axis. 


-M 

XX 

ydA  =     —  Q 

I       ap 

X 


But 


Hence: 


where  Q       =  noaont  of  area  from  a  to  p  about  x  axis, 
ap 

dM 

~     =  V  (!•©•»  shear  force) 
oz 


dP 
dZ 


I  X  dZ  \    I     / 


X  X 

Note:     The  second  tern  in  right  side  of  abore  equation  equals  zero  for 
non-tajwred  vlng. 
Therefore : 


ap 

q     +  * 


(1) 


Points  a  and  a'  as  ehovn  in  Figure  2  are  points  adjacent  to  either  side 

of  veb-slcin  Junction.  Figure  k   shows  the  spanvise  forces  acting  on  section 


aa'p'. 


P-niEdz 

<i2 


Figure  5 


From  Figure  6: 
^a  *  ^a-  -^  V 


dP  _   ap' 
dZ  "  I 


X 


(2) 


Considering  Bpanwia©  forces  acting  on  section  a*p'*  (Fig.  2  euid     ),  the 
equill'brlia  equation  glTos: 


•        3       Idzl  I 

\      /a'p"  X 


or 


S  =  \'  ^ 


^a- 


(3) 


Figure  5 

Considering  spanvlse  forces  acting  upon  section  p"  g'  (Tig.  2  and  5)^ 
the  equilibriuB  equation  glTes: 


h'*^~-  K 


from  which 

^g*     ^      I 


P  6 


P   g. 


Figure  6 

cmd  froM  equation  5*     q   ,  "  <1  ,  "••  ~Z' 

g  a'         I 


-  T  ^«  I  T 


(B) 


X  Z  X 

Considering  spanvlse  forces  acting  upon  section  gp"'  (Fig.  2  and  6 ),   the 
equillbrluH  equation  gires: 


dp 


^v::-. 


Sr2  "  ^«"  -  ^«'  '  - 


dP 
dZ 


gp 


"  I 


or  from  equation  B, 


v^L  ^'^L". 


Figure? 
Considering  spanvlse  forces  acting  upon  section  g'p*^  (f'ig*  2  and8  )^ 
the  e(iulllbrluii  eqiiatlon  girea: 


or 


^       gp  ^ 


(D) 


Figure  8 
Considering  apanvise  forces  acting  upon  section  p'^  d  (Fig.  2  and  9  ), 
the  equilibrium  equation  glres: 


S 


or 


-  ^ 


dP^ 
d-ldz,^p.^ 


VQ...^ 


!^ 


%»  =  ^d  "  T 


i^) 


Figure  9 


From  eg.uatlon8  D  and  k: 


n 


^6"        MI  I  *^d  I 

XX  X 


gfed 


And  equation  C  'becoiMs: 


S»2  "  '  '^d         I  \»         I  I 


J.      gP       ^      agred 
S#2       ^a»       MI  I 

X  X 


(5) 


Considarlng  spenwiso  forces  acting  upon  section  dd'p"'  (Fig.  2  and  10)^ 
the  equilibrium  equation  glYes: 


V  *  \  -  ^d' 


dP 
dZ 


dp 


n  I 


(E) 


Figure  10 
frotx  equations  ^  and  E: 


VQ^  VQ^  VQ^ 

M'       ^a'  I  I  I 

XX  X 


or 


VQ 


\'  "  ^a'  I 


(F) 


Considering  spanvlse  forces  acting  upon  section  d*p^  (Fig.  2  andu), 
tbe  equlll'brlum  equation  glyes: 


8 


%  -  ^d- 


dp 

dZ 


!«▼ 


d'p 


and  from  eq.uation  F 


VQ 


agfad      d^p^       ,^. 


Figur»  11 
Thus  eqiiatlona  one  throu^  six  reprosant  ths  spanwlse  equilibrium  equa- 
tions. These  equations  InrolTe  the  unknowns  q^,  q^,  q^  ,  <i^,  a,   q  , 

q  ,  q   and  q_,.  Four  additional  equations  are  now  obtained  by  conslder- 
a   a      d 

Ing  the  rate  of  twist  of  each  cell  and  a  moment  equation.  The  moment 
equation  vlU  Introduce  a  new  unknown  {i.e.,   the  shear  center  location). 


•    1' 


0  =  controld 
B.C.  =  shear  center 


Figure  12 
Referring  to  Figure  12  ^tnd  taking  moments  about  point  h: 


(G)  V{r  +  1    )+2  /   q  dA  -2 
B.C.    /  ^  n  If 

abc      _/  ag 


"b^cI-^/,  S'^BX-'f  «c"o2-V  .^^'^'^J  "  ° 
y  gf  ed      y  dc      y  gd 


Substituting  Into  equation  G  ralues  of  q^,  q^,  q^,  q  ,  q^  In  terms  of 


V  \*'   ^^  "^d* 


H.      V(r+x  )+2 

8.C. 


ap 
q   +  - — * 

V      ^        I 
abc  X 


f 

vq".  ,. 

r 

dA  -2 

4    .*-~^ 

dA      -2 

N 

^a'        I 

cl 

Ms 

X 

J 

gl'ed 

'd       I 


"r: 


-2 


Jdc\ 


q    .-.  -^5^  ^  -^^  I  dA  .-2 


c2 


J 


g«i\ 


M       I        .       ^a'       I  I  , 

Z  X  X  ' 


dA 


c5 


DlTldlng  aquation  H  by  V  and  recalling  q  ^   q   ,,  and  q,  are  consteuita: 

a       a  d 


2q  2q    ,  2q 

a  a  d 

I.      (r  +  X         )   =  -  A  + A     +  A„  + 

^  8.C.  V     n       V         cl       V      Rl 


2q  2q 

-^A      --^A      -^ 
V         c2       V       c3     I 


'J 


abc 


X   , 

Q     dA 
ap     n 


^7 


Q^        dA      -^ 
a'p"     cl     I 
ag       '  X 


( 


gf«d     ^  xy 


do 


Q  dA       +  A 

agfedg     c2         V         c3 


.c<'P^-^t)/«^-"^*t, 


I  .  2 

Q     d.<l     +  ~ 

gd  xy 


Sd 


Q         dA 


But  A,   +.A^  +  A^  =  A     (area  of  center  cell) 
cl         c2         c5        c 


And 


A^,  -  A   =  A  (area  of  rear  cell) 
Rl    c5    H 

Hence j  equation  I  becomos: 


2q. 


I.      (r  -»-  X         )   = 

8.C. 


2<1. .  2q 


^a'  d  2       X  r  T       2       X 


2_ 

I 


X  2 

ap     N     I 


.  I  abc  X  i  ag  x  I  gf  ed  x  J±c 

The  equations  of  the  rate  of  twist  for  each  cell  give  three  additional 
equations.     Note,   the  rate  of  twist  for  each  cell  Is  zero  since  the  load 
Is  through  the  shear  center. 


Front  Cell: 


J.    2a 


•\ii| 


abca 


J 


ds 
abc  / 


ac 


10 


K.   2GA 


Center  Cell: 
dQ  I 


c  dZ 


acdg   ^ 


da 
ac        J 


cd 


ds 
'c  t 


J 


g<i 


ds 


ag 


s 


ds 


=  0 


Rear  Cell: 

I  dQ 


dZ 


defg 


J 


gfed      J 


Sd 


V 


ds 


=  0 


Substituting  into  equations  J,  K,  and  L,   the  Talues  of  a  ,  q^^,  %.-t*   ^  ' 
a  ,  and  a  In  terms  of  q.  ,  q  ,,  and  q  , 


Fixsnt  Cell: 


r 


7 


VQ 
ap 
q  + ^ 

abc  ^      X  -< 


ds 
t 


J 


ac 


VQ 


-^- 


\'"  I 


3C  -» 

ap' 


X  -J 


d3 
t 


=  0 


II. 


<1   ^ 
a 


t  "*"  V 


J abca  J 


ds   1_ 
t  "^  I 
ac        J 


X   ds 

ap  t 


abc 


X   da   ^ 
%         --  =  0 
ap'  t 


7 


ac 


Center  Cell: 

r 


y 


ac 


X  — , 

VQ  , 
ap' 

q  -  q  ,  +  -r"^ 

a    a'    I 

X  J 


da 
t 


J 


cd 


a'     I        I 

X  I 


ds 

t 


"^^^•^         VQ^    VQ^„,' 


q  +  — ~ —  +  q  . 

^d    I      ^a'    I      I 


76d 
or 


X   -« 


ds 

t 


\'    I 


^   ^  P 


yag 


X   -J 


ds 


=  0 
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III.  -  zr 


ac 


ds 
t 


\' 


J 


ds   ^ 

^   t"  "*■  V  , 
acdga        J 


da    agfedg 


gd 


y 


da    agfed 
cd        I  y 


Sd 


ds 
t 


y  ac 


ap'  t 


X    de   1_ 
\'p'^  t~  ■  I 


V 


cd 


I    da   1_ 

OT'"  t  ■  I 


y 


8d 


^3c    da 
Q     —  =  0 
a'p"  t 


J  atj 


Boar  Cell: 


J 


gfed  '- 


VQ- 


X    -1 


IT 


^ 


X   -J 


de 
t 


J 


gdL 


VQ       VQ 
^       S^od  ^  __ag  ^ 


I 

X 


ds 


or 


IV. 


ds 

t 


^a' 


->'gfodg 


ds 
t 


agfed 


7  gd 


da 
t 


.>gd 


da 


*dp'^  t 


V 


gfed 


X  da 

^ag  t 


=  0 


J 


gd 


Thus  equations  I,  II,  III,  and  IV  represent  four  e(iuatlons  with  the 

four  unloiovna  q.  f  q.   ,,  q,.,   and  x    .  Once  these  eg.ua tlons  the  desired 
a   a    d       s.c. 

shear  flovs  can  he  obtained  by  uae  of  eqxMitlons  1  through  6.  These 
equations  then,  apply  to  the  solution  of  the  shear  flow  problem  InrolT- 
Ing  a  structure  of  three  cells,  skin  fully  ef fee tire  In  bending,  and  the 
load  acting  through  the  shear  center. 

To  the  shear  flows  prerlously  calculated  for  the  load  thix)ugh  the 
shear  center,  those  shear  flows  due  to  a  torsional  load  must  be  added. 


Figure  13 


s.c 


s.c. 


s.c. 


=  Shear  Center. 

=  Applied  Torsional 

Mooent  about  Shear 

Center. 
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From  the  prorlotiB  apanwlse  oqulllbrium  equations  and  the  fact  that 

dp 

-*  =  0  for  pure  torsion  the  following  equations  result; 

dZ 


Taking  noaaente  ahout  shear  center  (Fig.   1^)   glres: 


(m)       M  +  2 

8.C. 


dA.    +  2 


J 


ahc 


%^1 


-  2 


J 


J 


cd 


a  A  ^  -  2 

c  c5 


ac  / 

gfed 


.eu  _  -  2 


''v2*olt 


Suhstituting  into  equation  (m)   the  ralues  of  q^,  q^  ,  q^,   q^,   q  ,   and 

q^  in  ter«B  of  q  ,   q   ,>  and  q     glTes: 
x\  a       a  d 


M  +2 

s.c. 


abc 


q  dA,   +  2 
a     1 


(-q.q     )dA^      -  2 


ac 


'a     a'       c2 


J 


ag 


a'     c5 


-2 


J 


gd 


^^a'-^d^^cl.  -  ^ 


y 


cd 


a'     c5 


y 


gfed 


^d^6  =  ° 


But  A,    -  A^  =  A       (Area  of  nose  cell) 
12         W 

A^+A       +A,+A       =A       (Area  of  center  cell) 
c2         c5        c4        c5         c 

A^  -  A,    =  A       (Area  of  rear  cell) 
6        1+        R 

Hence : 


(V)      M  =  2q    ,A     +  2q  A     -  2q  A^ 

^  s.c.  a'  c  d  E  a  W 
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Considering  now  the  rate  of  twist  In  each  cell  axid  assunlng  that  the 
rate  of  tvlst  for  each  ceH  Is  the  aaao   glrea: 


Noae  Cell: 

dOJ 


(n)    2QA, 


N 


dz/ 


abca 


ahc 


ds 


J 


ac 


ds 


Center  Cell: 


(P)   2a^ 


dQJ 


\      /acdga 


; 


ds 

Sri  7- 
ac  / 


cd 


ds 


J 


gA 


ds 
^w2  t 


J 


ag 


ds 


Rear  Cell: 


'gdefg 


) 


ds 
gfed       J 


gd 


ds 

\f2.   t 


Substituting  Into  equations  (m) ,  (p),  and  (q),  the  values  of  a  ,  a  ,  q,^, 
q^  ,  q^,  and  q  In  terns  of  q  ,  q  ,  and  q  glres: 


Rose  Cell: 


(VI)    2GA. 


'd0\ 
N  Idz) 


=  1. 


abca 


abca 


ds 


ac 


ds 
t 


Center  Cell: 


\   /acdga 


Rear  Cell: 


J 


acdga 


ds 

r  -^a 


ds 

—  +  q 

t    ^d 


ac 


gd 


ds 

t 


/  dQ 

(vni)  2^\- 


=  -  <1. 


\   /gdefg 


/ 


gfedg 


ds 


/ 


gd 


ds 
t 


Thus,  equations  (V),  (VI),  (VII),  and  (VIII),  represent  four  equations 

/d©\ 
and  four  unknowns,  q  ^  q  ,*  q^>  aiid  — }  .  Once  these  equations  are 

a   a    d      ^dZ  \ 

solred  the  shear  flovs  due  to  a  torsional  load  about  the  shear  center 


W 


nay  be  dotormlnod  by  us©  of  ©qxaatlons  (7)  through  (12). 

Suanarlzlng  the  equations  for  the  deterailnatlon  of  the  shear 
flovs  In  a  three  cell  bean  vlth  skin  fully  effectlTe  In  bending  glTes: 

A.  Equations  for  Vertical  Load  Through  Shear  Cen'tor: 

.     ^^•^W   ^^*'^c   Vr   2   X     r^     ^  -r 


2  ^x  .  2 

+  —  Q   A  -  — 

I  ^gd  c2  I 

X     ^  X 


X  ,  2 

Q  dA^  +  -- 

ap  If   I 


J 


abc 


J 


q'  dA   -^ 
ap"  cl   I 

ag  j  gf od 


%.V^i 


2  r 


Q"^  -  dA 
^d'p^   c2 


;dc 


(II)  ~ 


yabca 
<1-  '' 


(III)  -  f 


da   V 
t  "*"  V 


ds   V 
t  "  V 


de   1 
ac      z 


Q^      ^      I. 

V    *P   *   '  I 

abc  I 


/ 


ac 


«3C   de 
ap'  t 


J 


ac 


y acdga 


t  "*"  V 


de   ^x 

—  -  Q 

t     a^edg 


0<i 


'  da 
t 
/cd 


agfed 


ds   1_ 
t  "*"  I 


8d 


X   ds   1^ 
ap'  t  "  I 


ac 


\'P'V   t 


Vcd 


X    ds   1_ 
%•"  t  ■  I 


ffi 


«3C    ds   _ 
0     —  =  0 
a'p"  t 


ag 


(IV)  !d 


ds 


Q  y  Q 

a'f   ds    agfed 


t    V 

ygfedg      y 


ed 


ds   1 


6d 


X     dfl 

%•▼  t 


y 


gfed 


X  ds 
ag  t 


yad 
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(1)  <i.=  'i.*v«^^ 


<.. 


(2)  S,x=-'a-V*-f^ 


™.^ 

^..^" 


^5^  S  =  ^a>"i" 


VQ 


('^^  %.  =  ^d--i 


d£ 


nr 


,^v  agfod  ^   gp' 

Z  X 

(6)   q^  =  q^,  -.  -- f-*  -^  -J^ 

X  X 


B.  Sguatlone  for  Torsional  Load: 

(V)   M     =  2q  A  +  2q^A^  -  2q  A^ 
•.c.     a'  c     d  R     a  H 


(VI)  2GA, 


dC 

dzL   "  *^a 
\   ' abca 


r 


da 

—  +  q 
t    ^a' 


J 


abca 


da 
t 


J 


ac 


(VII)  2GA 


de 


acdga 


-'\, 


da 
T-^a 


(VIII)  2GA^ 


dO 
dZ 


acdga 


J 


d8 

—  -♦•  q 
t    M 

ac       J 


(      ds 
t 

Si 


=  -<1. 


gdefg 


ds 
T^^a' 


t 


J 


gfadg 


J 


gd 


(T)    «n  =  «« 

(8)    ^^  =  -.^. 

•».' 

(9)    ^  =  v 

(10)      ^   =   q. 
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(11)  S2  =  <1».  -1, 


(12)  <l^  =  <l.. 


<^"   111     =  i      -  (I 

\   fabca      facdga      /gdefg 


R»f erring  to  the  above  euauirj  of  equations,  the  resulting  shear 
flovs  for  pure  torsion  are  of  constant  Talue  betveea  any  tvo  success Ito 
structural  Junctleas^  vhlle  the  shear  fletrs  due  te  bending  rarj  betveea 
such  Junction  x>ointe.  The  ■stIiw  and  ■1nl»u»  ralues  of  these  shear 
flows  can  of  course  be  obtained  by  differentiating  the  sbsco*  flev  ex- 
pression vlth  respect  to  the  paraaeter  defining  the  aaount  of  area  of 
the  appropriate  section. 

The  procedure  oHplojed  in  the  aaaljsis  of  the  subject  three  cell 

beam  is  applicable  to  a  *n'  cell  beaa  (for  n  >  2)  vhere^ 

^  -  2  =   number  of  uakaovns  (i.e.,  shear  flovs  and  shear  center 
location) . 

n  -   Buaiber  of  equations  arailable  froM  the  rate  of  twist  of 
each  cell. 

1  ~  asBsnt  equation. 


Thus: 

(^  -  2)  -  (n  -  1)  =  2n  -  3  =  nuaber  of  points  at  which  the  shear 

Talue  is  assuMsd. 

FreTleusly  it  was  shova  that  when  considering  the  principal  axes 

'ap     X 

If  the  axes  are  not  principal  axes  and  in  addition  the  load  is  acting 
la  anj  giren  direction  then 
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p  = 


MI   -MI 
3C  xy    J  T 

I  I  -  I- 


X  + 


M  I 


yap 


X  7 


xy 


I  I 
X  y 


M  I 
X  y 

xy 


dA 


MI   -MI   ^ 
X  xy    y  X 

I  I  -  I^ 
I  y    xy 


y  ap 


MI   -MI 
X  y   xy 


7  ap 


ydA 


MI   -MI 

I  I^^ 
X  y   xy 


MI   -MI 

ap   I  I  -  I^ 


P  = 


I  I 
X  y 


X  y 


xy 


ap 


xy 


[". 


(I  Q^ 
xy  ap 


I  Q^  )  +  M  (I  Q^  -  I  Q' 
y  ap    y  xy  ap    x 


apj 


^   '  dZ    II 


X  y   xy 


V  (I  </ 
[^y  xy  ap 


-  I_Q^  )  +  V  (I  _Q^ 


X  dZ 


/l  Q^  -IQ^\ 

II  -  i^ 


X  y 


xy 


y  ap 


♦  M 


X  xy  ap 


dZ 


I  Q 


\ 


I  I 
X  y 


IQ^  ) 
X  ap 


X  ap 


xy  li 


Thus  the  algetra  now  InrolTod  Is  rery  cuiibersiaii,  and  a  practical 
analysis  would  require  eoae  assuaptlona  or  approxlBatlons  wherein 
this  difficulty  could  "be  orercone.  For  exaaple,  in  the  procedure 
employed  herein  it  vas  aseuaed  the  wing  or  beam  had  no  taper,  hence, 
the  last  two  terae  of  equation  (IX)  each  equals  zero.  Also  principal 
axes  were  used,  thus  I   equals  zero,  and  only  a  rertical  load  (i.e., 
perpendicular  to  the  x  axis)  was  considered.  These  qualifications 
reduce  equation  (IX)  to 

dPi   _   ap 


dZ 


\ 


ap 


Id 


Al69  vorthj  of  note,  Is  that  the  shear  center  location  as  derlTed 
herein  Is  different  than  that  location  as  obtained  from  the  analysis 
where  the  skin  Is  considered  either  non-ef foe tire  or  partially  effec- 
tlTe  in  "bending. 
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PART  II 


STRESS  ANALYSIS  AND  STATIC  TEST  OF  A  MODEL  WING 
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PART  II 
STRESS  AHALYSIS  AID  STATIC  TEST  OF  A  MODEL  WHIG 

This  part  of  the  report  Includes  the  stress  analjsls  sjid  static 
test  of  a  three  cell  aedel  vlng.  The  vlng  vas  constructed  from 
.O72-2UST  aluBlnua  sheet,  vlth  six  foot  span,  tvelre  Inch  chord,  and 
circular  arc  alrfell  shape.  As  far  as  the  stress  analysis  Is  con- 
cerned, those  equations  deteloped  In  Part  I  could  he  applied  directly, 
hoverer,  due  to  the  STsnetrj  of  the  vlng,  the  shear  flew  ec^uatlons 
are  soaevhat  slapler  than  those  of  Bart  I,  hence,  a  separate  stress 
analysis  for  the  vlng  using  the  prerleus  procedure  Is  Included.  The 
Initial  sections  of  Part  II  are  deroted  to  this  stress  analysis  and 
the  actual  static  test  follovs. 

Geoaetrlc  Properties  of  Model  Wing 


Due  to  Bjwmetrj   aheut  chord  line 
only  upper  half  of  vlng  need  he 
considered. 


^  +  x^  =  R^ 


(y+a) 

a+i  =  o,y  =  ,6k 

.61f4a  =  R 


t  =  .072" 
R  =  28.M*8" 
c  ^  12" 

a  =  R  -  .6U  =  27.805"  *  -  27.805" 
a  _  27.805 
R  "  2Q.hh^ 


sin  7 


=  .9775 
7  =  1.5585  rad. 


R^  =  809.118  sa.ln. 
a^  =  775.118  sa.ln, 
7  -  1.3585  rad. 
sin  7  *   .9775 
aR  =  790.9152 


Fig.l 


Itomut  of  Inertia  (I     ) : 


/  R  sin  0  =  y  +  ft 

Tr/2 
^i  Skin  '  ^    I  ^^  ^^'^  ^  '  *^^  ^^  ^ 


'7 


0 


2     ■    4 


=    liRt  I  -7^    -  ~    slB  2<^   +  2«R   CO8    0   +   R^^ 

1.5585 


772 


=  8.192  {UOI^.559  [1.570796  -  1.5585]  -  202.279  [0  -  sin  2.7170] 
+  790. 9152  X  2  [0  .  coe  1.5585]  +  775.118  X  .212296  ] 

=  8.192  {kOk.^'^k  X  .212296  +  202.279  ^  .Ull9'*9  +  158.18264 X 
[-.21070]!  +  775.118  X  .212296} 

=  8.192  [85.885  +  85.52865  -  555.29082  +  l6lf.l298} 

=  8.192   A. 05261  =  .1+509895  in.* 

I      V4     =  -^509895  in.* 
I  skin 

I  ,  =  2  X  —^  =  '^^   ^  •^''  =  .012  X  .875156  ^  .0105016 
wets      12        6 


flanks 


]  -  '  [■'  ^•'^''  *  -^  '  •'^^  "  •'^^5'] 


=  .0097648 


r  ^  ,  =  .MHU9  In.* 
total 
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GooBotric  Propertloe  of  Model  Wing: 


Mote:  The  rlret  line  for  the  spar 
Is  k   inches  from  trailing  edge.  The 
vel)  location  Is  3*^  Inches  from 
trailing  edge.  For  purposes  of  analy- 
sis the  spar  loading  Is  assumed  to 
transfer  to  the  skin  at  some  point 
between  the  rlTet  line  and  veb  lin« 
(«.g. >  5.79"  from  t.e.). 


.  i,/       2.21    2.21     ^,,     ^06     6 
8in(^  =  -^-   =  r,Q  W^   =  -0777    sin  P  =  -  =  - 


R   28.U45 
^  =  .2126  rad. 


SO  =  R  cosi//  =  28.M*5       .99698  =  28.5591 
aY  =  St  =   .55^1" 


COB  If/  =   .99698 


Rd0  =  R((/'  -  0)  =  28.UU5   .0778  =  2.215" 


Rd0  =  rO  -4^)  =  28.1^1^5  (.2116  -  .0778)  =  3.802" 


cos^  =  .9775 


/ 


ds  »  2  I   ds  =  7.604" 
ate        J ab 


ds  =  2 


jsa' 


ds  =  lv.l^50" 


J 


La' 


; 


ds  =  2  /     dy  =  1.10" 
ac       Jo 


1      Area  "aOL' 

7  A  :  rs 

k     c  R 


£ 

2TT 


or   Area  'aOL 


.^OL.  =  .°n8iJ09,ll8  ^  j^  ,^,59 
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Area  "aos"  =  J  X  28.3591  ^  2.21  =  31.3368 

^  A    =  2.21  X  .55  +  51.1^759  -  31.3568  =  1  3552  sq.ln. 
c 

A     =  5.i*128  eq.   In. 
c 

4  A     =  Area  ""boL"  -  -^  A     -   "tot" 
'^     n  c 

Area   '^oL*         P             .          ,^  ,„       .2126  X  809.II8       q^  ,,,,^, 
-* rr =  —  Area  "boL"  = ^ =  86.  OO92I+ 

Area  "bot"  =  Jx6xa=3x  27.8O5  =  83.U15 

J  A     =  86.00921+  -  1.3532  -  83.1*18  =  1.2lf07 

A     =  2.U8llv  sq.   In. 
n 

Calculation  of  Mo«ent  of  Areas:     (Refer  to  Fig.   6) 

9  r© 


J 


CL    „  ds;       *^  ..  =    /       y  Et  dO  =    /       Rt  (R  cos  6  -  a)  d  0 


La»       *^  '         /O  yO 


.1 
'La* 


=  Rt  [r  ein  Q  -  ao) 
^^=.0778 

R^t  [r  8 In  0  -  aoj  dO 


\p"  ^  = 


h 


-^   .0778 
=--  R^t  <  -  R  cos  0  -  —  r 

=  809.118  X   .072  ( -28. Ii45(. 99698  -  1) 
27.805       .0778^1 

=  58.2561*9  (.0859039  -    .O841U65] 
=   .10237995 
2U 


f^ 


Q       ds:       Q 
ap  ap 


^Q 


Rt  ydO  = 


yab 


-Jv 


Rt  (R  cos  0  -  a)   d9 


h 


=.  Rt  JB(8ln  9  -  Bln(f')   +  a((^-     8)] 


Q      ds  = 
ap 


'ab 


R^t  [R(sin  0  -  sln^)   +  a  (y/-  0)]] 


dO 


V 


,P 


rac 


v^     ,  de: 
ap' 


2     r                                                           aO^  1 
=  R  t     -  R  cos  Q  -  RO  sin       +  a       ©  -  

=  R^t  f-  R(coB  ^    -  cos^)    -  R  O  -y)   sinjt/  +  a^P 

a£f       alfl 
2     "       2  J 

=    56.25649  (-  R(-. 019^)  -  R(.  1548) X. 0777 

*  a  X  .^8  X  .2126  -  i^^  .  2LS^] 

=  58.25649  {.55^1086  -  .29795179  +  .45990248 

-  .0841496  -  .6285551] 
=  58.25649  X  .00557659 
=  .208559579 
h=.554l 


«ap-^ 


yt  dy  = 


2     2    2  ^     "  '' 


'ac 


t   .  ds 
ap' 


h=.554l 


=  2 


=  t 


['•  -  ^1 


\  (•>=  -  y") 


2h^t 


dy  =  t  Ih^y  -  y^/5 


0 


=  2JL.^2ip.^2I2  .  .00816592 
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ac 


Q^   da 
a*L 


<•. 


yRtdO  =  /   Rt(R  cos  0  -  a)  dO 


[■ 


=  Rt  R  sin  ©  -  aO 


-1.0778 

Jo 


=  2.04804  [28.445  /  .0777  -  27.805  X   .07783 
=  2.04804  [2.2101765  -  2.165229] 
=  .096147456 


ac 


Q  ,  de  =  .09614745  X  1.10 

a  L 


=  .10576219 

G©o«©tric  Properties  of  Wing 
(Strain  Gege  Locations) 


,!^ 


76.^68- 


/V,968 


i. 


, US 


,  :A 


bOH7  I 


I 


':^ 


(>.o3l 


\\ )[_ 


r^.  .^ 


jvy 


-i  1 


;6,<?07- 


26. Z3 


/v3a^ 


Lo<3d  Point 


Fig.  3  PlanforH  of  ving  showing  strain  gauge  locations : 
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Referring  to  Fig.    J,   the  strain  gauge  locations  are  those  ac- 
tually used  in  the  static  test.     They  are  Introduced  here  In  order 
to  analytically  determine  the  shear  stresses  at  these  locations. 


Stress  Analysis  for  Vertical  Load  Through  Shear  Center 


L   q-O 


p^f^n   ^" 


c  c' 


•K  <^»0  X 

S'^      I 


P.f^^l 


(1) 


VQ 
a  p 


Kv  ^p" 


'   *.•  =  T 


■s:  ^.-JiE: 


<'. 


(2) 


S,  *  1.  -  \' 


ap' 


^a         I 

X 

+  -- —       (3) 


YQ 


\  -    \  =      I 


SIL 


(^) 


rig.  4 


Twist  of  KoB9  Cell: 


ahc 


ds 

^nT 


J 


ac 


ds 


=  0 


(5) 
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(<1  +  -^)   48  - 


'a    I 


J 


4 


abc 


_/ac 


-Q   + 

^a    I 


ap*     a'L 


ds  «  0 


X  -^ 


de  +  q 


abc 


TQ 


J 


de  + 


ac     / abc 


±2 


ds  - 


J 


ac 


ap'     a'L 

I    "*"   I 
I-  X       X  J 


ds  =  0 


2a. 


de  +  q 


y  ab 


da  + 


J 


ac 


2V 


y  ab 


Q^  <ie  -  —  I   (Q^  ,  +  Q  ,J  de 
ap     I      ap'    "'T 


a'L' 


ac 


2V  V 

2q  (5.902)  +  q  (1.10)  +  -"(.2083596)  -  --(.OO81659  +  .IO576219)  =  0 
a         a        I  I 

X  X 

8.70if  q    +  7-  r  .*H67192  -  .1159281]    =  0 
a       1  -* 


8.70lv  q    +  .5027911  7- 

=  0 

X 

q^  =  -.05^79^ 

X 

q     =   -.0788  V 

a'L                 X 
\'         I                       a'L  " 

Rt  (R  C®B  9 

-  a)  dQ  =  Rt  [r  sin  y/  -  a  If/  ] 


Jo 


2.0it80U  [28.UU5  X  .0777  -  27.805  X  .0778] 
2.0480U  [2.210176  -  2. 165229 J 


=   .0469^75 


.Oif69^75 
\'  '  .U4U9 

q  ,  =  .1065  V 
a' 


28 


strain  G«ige  C. 


Figure  5 


^"^  ^         R  28.MV5      -^^^ 


a  =  .OUfilt  rad. 


of  „  =  Bt  [r  Bin  0  -  a©] 

of     =  2.0l«£0»t  [28.1^45   y  .01^3  -  27.805  X  .0481*] 

=  2,0kQ()k  [1.37369  -  1.5^37623 
=  2.01*601^  ^  .05015 


'] 


I 
gage  c  X 

Tc=  .15977  V 


7"  ,  —  V  =  1.9^12  Y 

c        t 


=   .06171 


From  eq.t.    {h)  pa. 


VQ' 


©r 


%  - 1.  = 


aB 


'aB 


f 


Bin 


\  R 


2Q.kk^ 
p  =  .1679 


.1669 


J 


Rt  [r  cos  0  -  aj  dQ 


0788V  +  .28O7T 


V 


q^  =   .2019T  y    =  2.80»tY 
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=  Rt  [R(8ln  p  -   sln(//)   -  ap  +  a^ 
=  2.0l^0U  [28.Mf5(.l679  -   .0777)] 

-  a(.l679  -   .07783 
«  2. 01^804  [2.565759  -  2.5052505] 
=  2.01+801*  X  .06051  =  .1259269 


BexsLce^   the  shear  stresses  ae  aoalytically  deterKlned  at  strain 

gauge  locations  B  and  C^  due  to  a  Tertical  load  through  the  shear 

center  of  the  vlng  are: 

At  lecatlon  B:     T    =  2.80i*  V     p.s.i. 
'B 

AtlocatienC:      T    =  1-9^1  T     p.e.i. 
'c 

where  V  =  applied  load- 


Stress  Analysis  for  MoBont  about  Sheeir  Center 


\.   -%-^^     \-%       (*) 


a' 


<i^*\-\,  =  ^  s,-\^-\     (^) 


-  a        (c) 


Due  to  Bjm^trj  q     "  ~  "Id         ^^^ 


Figure  7 
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Taking  MoMonts  about  Point  0: 


a  dA  .  +  2 


cl 


La< 


±k  .   +  2 


Si^nl 


J 


abc 


\ 


c2 


J  ac 


q  dA 

M   c3   2 


c'k 


\' 


dA  ,  +  q 
cl    a 


J 


La' 


"nl  *  \' 


1  ) 

a 


abc 


J 


c2   ^a* 


dA 


c5  "  U 


ac 


7 


a'k 


\' 


dA   + 
cl 


;La' 


y 


dA  ^  + 
c2 


ac 


dA 


c3 


c'k 


+  <1. 


dA 


nl 


abc 


dA 


C2 


ac 


T 
k 


(5)   <1,A  +q  A  =7    whor©  A  =  4  area  of  center  cell 
a'  c    a  n   4  c   ^ 

A  =  area  of  nose  cell, 
n 


Twlat  of  Nose  Cell: 


2GA 


d©  _ 
dZ  " 

ds 

ds 
S.T 

yabc                 ^ac 

'\ 

y  abc 

t 

f 

ds  + 

de     - 

) 

abc           / 

ac 

d©  ^  _ 
dZ  '  i 

1 
>GA 

1 
a 

ds  ■^  q 

da 
t 


/  ac 

q  f 
^a' 


ds 


abc 


/  ac 


4b  -  1^. 


ac 


ds 


ac 
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Twlet  of  Center  Cell: 


fdOJ 

f 

2GA 

= 

c 

dZ 

/i 

da 

ST* 


d8 


J 


ac 


ds 

ST* 


J 


cd 


da 
SrT 


J 


dg 


2  ' 
.t 

7 


q^ds  .  I  I    q^  ds 


ga 


ao 


2 
t 


\'^^  *  t 


y 


ga 


(q  ,  -  q  )  d8 

a     a 


y  ac 


d© 
dZ 


2 
t\. 


GA  t 
c 


ds  + 


_7  ga    y  ac 


ds 


^\^ 


ds 


y  ac 


*«■ 


d.  +  1^. 


ds  -  q 


J  ac 


ds 


J 


ac 


dZ 


nose 
cell 


dO 
dZ 


center 
cell 


-7-  fq  (T.60lf  +  1.10)    -  q    ,X    1.10]     =  7^  fq   ,(^.**50  +  1.10)    -  q    X    l.iol 
fcA     l_a  a  -JAia  aj 

n  c 

"5  U128 

^x   2.U814  (^-7^  ^a  -  ^-^^  ^a'^   =  5-550  q^.    -  1.10  q^ 

^.7^655  q     -   .599865  q   ,  =  5.550  q   ,   -  1.10  q 
a  a  a'  a 

5.84655  q    =  6.12986  q  , 
a  a 


(6)  q     =  1.01|64  q   , 

a  a 


7roH  eqt.   (5) 


5.4128 


5.50789  q^,  ^l 


q   ,   +  2.U814X  1.0U8l+q   ,  =  - 

q   .  =   .OU71  T  T  .  *  -6542  T 

a  •  '  a  • 
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a  a 

q     =   .OU938  T  T     =   .6858  T 

•  a 

It 
Hence,   the  shear  stresses  as  analytically  deterBilned  at  strain 
gauge  locations  B  and  C,  due  to  a  BOBent  al>out  the  shear  center  of 
the  wing  are: 

At  location  B:   T     =   .6858  T     p.s.l. 
B 

At  location  C:    ^    =   .651^2  T     p.e.l. 
where  T  =  applied  aoaent. 
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SQUIPMEJT  AID  roOCKDURE 

The  gooBotry  and  aaterlal  of  the  Model  vlng  has  heen  discussed 
prerlouslj.  An  atteapt  vas  Bade  to  test  the  vlng  in  pure  bending 
(l.e.^  no  torsion)  and  In  pure  torsion.  Six  SR-k   strain  gages  vere 
■ounted  on  the  vlng  and  an  additional  one  Installed  upon  a  saall  sheet 
of  2UST  alumlnuB.  for  temperature  correction  purposes.  (See  Fig.  5  for 
strain  gage  locations.)  The  rig  by  vhlch  the  bending  test  vas  con- 
ducted Is  shovn  In  Fig.  9-  Referring  to  Fig.  9>  the  load  vas  applied 
hj  adjusting  the  tumbuckle  vhlch  appears  on  the  right.  The  aoount 
of  load  vas  determined  by  a  prerlously  calibrated  dynanometer  installed 
Just  belov  the  tumbuckle.  Strain  gage  readings  vere  obtained  at 
sereral  load  ralues,  Tarylng  from  zero  to  kyj   pounds.   (See  Tables  I 
and  II  at  end  of  this  report.) 

Fellovlng  the  bending  test  the  vlng  vas  tested  In  torsion.  The 
rig  used  here  Is  shovn  In  Figs.  10  and  11.  Referring  to  Fig.  10, 
the  load  vas  applied  by  a  hydraulic  Jack  and  amount  of  load  deter- 
mined by  use  of  the  dynamometer  used  In  the  bending  test.  Not  ap- 
parent from  Figs.  10  and  U  Is  the  balance  scheme  used,  vhereln  a 
cable  vas  attached  to  the  free  end  of  the  vlng  and  loaded  vlth  sand 
bags  such  that  the  vlng  vas  not  subjected  to  bending  due  to  Its  ovn 
velght.  Strain  gage  readings  vere  obtained  at  several  torsional 
load  ralues,  varying  from  zero  to  over  6000  in. lbs.  (See  Table  III 
at  the  end  of  this  report. ) 


3^ 


i 


Fig,9.   Bending  Test 


I 


Fig.lO.   Torsion  Test 
o^       


^ie.ll.  Torsion  Test 


RESULTS  AUD  DISCUSSIOI 

From  the  resulting  test  data  (l.e  ,  Tables  1,   II,  and  III)  itohr 
stress  circles  were  drawn  for  each  strain  gage,  corresponding  to 
V  =  400  pounds  in  the  bending  test  and  T  =  625O  inch  pounds  in  the 
torsion  tost.  This  was  done  somewhat  as  a  preliminary  check  (i.e., 
a  comparison  between  experimental  shear  stress  ralues  with  those 
calculated  previously).  Considerable  disagreement  in  results  was 
Indicated  at  this  point.   It  appeared  that  perhaps  the  bending  test 
actually  introduced  some  torsional  moment  upon  the  wing  and  similarly 
the  torsion  tost  Introduced  some  bending  upon  the  wing.  Baring  suf- 
ficient data,  it  was  possible  to  attempt  to  separate  the  torsional 
shear  stresses  from  the  bending  shear  stresses  in  both  the  bending 
and  torsion  tests.  This  was  done  and  the  experimental  results  were 
coiqwnsatod  for  accordingly.  HowoTor,  there  was  still  considerable 
disagreement  between  the  experimental  and  analytical  results.  The 
stress  circles  are  not  Included  in  this  report  since  their  enly  ap- 
parent Talue  was  to  show  considerable  disagreement  between  the  ex- 
perimental and  analytical  results- 

In  Tiew  of  this  disagreement  in  results,  the  next  logical  ap- 
proach appeared  to  be  an  attempt  at  determining  the  reliability  of 
the  experimental  strains.  Figure  12  is  a  plot  of  load  rs.  strain 
from  strain  gages  C,  B,  and  E.  From  this  plot  it  is  apparent  that 
there  was  no  OTor-all  linearity  of  stress-load  Tariation  as  pre- 
dicted by  Hookes  Law.  The  causes  leading  to  this  non-linearity  be- 
tween stress  and  load  might  be  attributed  to  the  inherent  random 
error  in  the  strain  determination,  instability  of  the  structure 
wherein  the  slippage  at  the  riveted  Joints  causes  portions  of  the 
structure  to  load  or  unload  in  some  unpredictable  fashion,  error 
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du©  to  local  affocta  of  the  rirots  (i.e.,  St.  Vonant  principle), 
Bcalo  factor  of  the  nodel  and  limitations  of  the  strain  gages. 

In  discussion  each  of  the  abore  sources  of  error,   the  results 
of  Figure  12  indicates  the  inherent  randoa  error  of  the  strain 
dotorMlnations  was  of  ainor  Importance,   since  the  irregularities 
in  the  strain  are  of  much  larger  magnitude  than  those  attrihutablo 
to  random  error,  where  It  Is  known  that  the  strain-meter  readir^gs 
as  such,  are  at  least  accurate  to  within  plus  or  minus  one-half 
mlcroinch.     Error  in  the  load  determination  may  he  ruled  out  since 
such  errors  would  need  he  in  order  of  magnitude  of  plus  fifty  pounds 
at  some  points  and  minus  fifty  pounds  at  other  points. 

Errors  due  to  scale  factor  of  the  model  and  limitations  of  strain 
gages  appear  as  possibly  large,  when  considering  that  the  magnitude 
of  shear  stresses  is  quite  low  and  the  shear  stress  values  so  critic- 
ally dependent  upon  the  precise  location  of  the  strain  gages.     For 
example,   in  bending  with  the  applied  load  equal  to  approximately 
four  hundred  pounds  the  shear  stress  in  the  skin  orer  the  center 
cell  of  the  wing  raries  from  zero  at  the  center  of  the  panel  to  ap- 
proximately eight  hundred  pounds  per  square  inch  at  the  web  skin 
Junction.     Thus  the  shear  stress  reuries  from  zero  to  eight  himdred 
pounds  per  square   inch  orer  a  distance  of  about  two  inches,   hence, 
the  requirement  for  precise  location  of  the  strain  gages.     In  this 
respect  It  should  be  noted  that  when  using  a  strain  rosette  the 
strain  readings  are  assumed  to  represent  those  strains  at  the  inter- 
section of  the  gage  axes  on  the  strain  gage   (Fig.   13). 

Figure  13  is  a  full  scale  diagraj»atlc  representation  of  the 
strain  gages  used  in  this  experiment.     Consider  any  gage  axis  (e.g., 
gage  axis  3)  and  suppose  that  wb«n  mounted  on  a  structure  the  strains 
In  the  3  direction  hare  considerable  rariation  oyer  the  gage  length. 

4C 


Fig.  15 
Th©  strain  gag©  reading  in  th©  3  direction  would  indicate,  approxl- 
nat©ly,  an  average  strain  orer  the  gage  length,  and  as  previously 
stated,  it  is  assuaied  that  the  resulting  strain  is  representative 
of  the  strain  at  the  gage  axes  intersection   However,  it  is  quit© 
possible  that  a  truer  location  of  th©  r©3ulting  strain  would  b© 
soaowhat  to  th©  right  of  th©  gag©  8ix©3  int©r8©ction  point,  ©specially 
if  th©  l©ad  in  wires  to  th©  gag©  have  a  portion  glued  to  th©  sur- 
face of  th©  structure.  From  close  observance  of  Figure  15^  it  ap- 
pears that  the  true  location  of  the  resulting  strain  aay  be  as  nuch 
as  a  quarter  of  an  inch  to  the  right  of  the  gage  axis  intersection 
point.   In  usual  applications  of  strain  gages  this  error  would  be 
n©g].lglble,  however,  as  previously  Mention,  in  th©  3ubJ©ct  t©st  th© 
strains  might  vary  considerably  over  the  distance  of  onle  one-quarter 
inch.  This  error,  if  present,  cannot  be  coapensated  for,  since  if 
th©  tru©  strain  locations  along  th©  ax©s  are  not  at  th©  intersection 
of  the  axes,  then  data  taken  from  the  gage  r©ading8  would  b©  U8©d  for 
constructing  on©  strain  circle  to  represent  three  different  points 
on  th©  structure,  which  is  moemingless.  But  such  error  can  be  dis- 
counted as  the  Bain  source  of  error  in  this  ©xperlment,  sine©  such 
an  ©rror  would  not  caus©  th©  abinipt  changes  in  slope  of  the  strain 
vs.  load  curves,  in  fact,  the  ^rror  would  be  linear  with  load.  Also, 
there  is  present  th©  fact  that  such  ©rror  could  not  account  for  the 
discr©panci©3  in  th©  torsion  test,  since  here  the  shear  flows  are 
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constant  over  any  integral  section  of  a  cell  ("by  Integral  section 
is  Meant  the  skin  as  opjweed  to  the  webs). 

The  error  due  to  the  local  effect  of  the  rivets  could  certainly 
"be  prominent  for  soiae  of  the  gages  "but  would  not  explain  the  discrep- 
ancies shown  "by  gage  B  which  was  located  at  the  middle  of  the  center 
cell  skin  and  at  maximum  distance  from  a  rivet  line. 

The  insta"bility  of  the  structure,  wherein  slippage  at  the  riveted 
Joints  causes  portions  of  the  structure  to  load  and  unload  in  some 
unpredlcta'b'le  fashion  remains  as  the  possibility  of  the  cause  for  the 
large  variations  in  stress-load  slope.  A  further  study  of  Figure  12 
makes  this  appear  as  a  reasonable  conclusion.  For  example,  note  the 
strains  at  the  loads  of  U25,  565,  and  25O  pounds.  Ideally,  these 
strain  values  would  form  a  straight  line  for  each  gage  axis,  but  if 
a  straight  line  is  drawn  between  the  25O  and  h2^   pounds  strain  values 
the  363  pounds  strain  value  does  not  lay  upon  this  straight  line  for 
emy  of  the  gage  axes.  The  5^5  iwunds  strain  values  do  though,  show 
a  uniform  trend  upon  considering  that  gages  G,  B,  and  C  have  the 
same  gage  axis  direction  orientation  on  the  wing,  and  the  3^3  pounds 
strain  values  all  lay  to  the  left  of  the  aforementioned  straight  line 
for  gage  axes  Bi,  Ci,  Gi  and  B2,  Cg,  Gz,   and  to  the  right  for  gage 
axes  B3,  C3,  cuid  G3.  Here  again  is  evidence  that  the  error  Is  not 
caused  by  either  random  error  or  error  in  say  any  one  set  of  gage 
readings,  but  rather  that  the  structure  was  behaving  as  the  data 
indicates,  such  behavior  being  Induced  by  slippage  of  the  riveted 
Joints. 

The  slippage  of  the  riveted  Joints  could  be  expected  to  not  be  - 
continuous  with  load,  but  intermittent,  such  that  portions  of  the 
structure  would  "lock"  over  certain  load  ranges,  slip  and  lock 
again.  Thus,  it  is  probable  that  over  some  loading  range,  during 

42 


vhich  the  structure  remained  in  a  locked  condition,  there  would  be 
the  expected  stress-load  linearity,  and  such  is  the  case  for  the  load 
range  of  128  to  250  pounds.  Hence,  it  appeared  that  for  som»  load 
rang©  the  test  data  might  "be  representative  of  the  structure  without 
slippage  of  Joints.  To  determine  this  range  the  plot  of  Figure  l4 
was  constructed.  The  analytical  values  of  strain  as  shown  there  are 
Included  in  Table  IV.  From  Figure  Ik,   it  appeared  that  the  most 
probable  load  range  in  which  the  experimental  strains  would  repre- 
sentative of  the  actual  wing  without  slippage  was  between  l66  and  2^0 
pounds.  Even  in  this  range  the  experimental  and  analytical  stress- 
load  slopes  are  not  in  agreement  but  this  can  be  explained  from  the 
fact  that  even  though  the  test  was  intended  for  pure  bending  there 
was  some  torsion  entering,  as  indicated  by  the  strains  from  gages 
0  and  F.  The  next  step  was  to  determine  if  the  amount  of  torsion 
entering  would  compensate  for  the  difference  between  experimental 
and  analytical  stress-load  slopes.  This  was  done  from  analysis  of 
Figure  15- 

From  Figure  15  note  the  experimental  load- strain  line  for  gage  C. 
This  was  corrected  for  torsion  by  subtracting  those  strains  due  to 
torsion  as  indicated  by  the  shearing  strains  from  gage  G.  The  re- 
sulting load-strain  lino  showed  close  agreement  with  the  analytically 
derived  line.  However,  the  same  corrections  when  applied  to  the 
strains  of  gages  B  and  D  do  not  give  agreement  between  the  experi- 
mental and  analytical  results,  hence,  the  indication  that  slippage 
occurred  over  this  load  range  also. 
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LOAD  vs  SPIEAH  STRAIN 
EXPERIJENTAL  <?:  AJ^TALYTICAL  RESULTS 
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COHCLUSIOHS  AMD  RBCO>«ffiHDATIOHS 

The  results  of  tba  static  test  of  the  vlng  are  of  little  Talue 
but  do  Indicate  soae  recornnendatlons  for  such  a  test.  Jor  such  a 
test  the  scale  of  the  nodel  should  "be  larger  In  order  that  the  "life 
size"  rivets  do  not  predominate  upon  the  stress  distribution  of  the 
model,  and  also  to  reduce  the  preclseness  requirement  of  the  strain 
locations,  as  previously  discussed.  Also  care  should  he  taken  to 
proyld,e  rigid  Joints,  so  that  the  slippage  effects  as  evidenced  hj 
this  experiment  would  be  minimized.  In  fact,  It  appears  that  a  dif- 
ferent aothod  of  Joining  the  skins  at  the  leading  and  trailing  edges 
would  be  preferred  to  riveting. 

Another  reconoendatlon  Is  to  have  the  spar  flanges  as  small  as 
minimum  rivet  edge  distance  would  permit.  This  conclusion  arises 
a  condition  whereby  the  load  transfer  between  spar  and  wing  skin 
does  not 'necessarily  take  place  at  the  rivet  line.  For  example, 
suppose  that  the  skin- spar  force  relation  Is  such  that  the  two  sur- 
faces are  being  pressed  upon  each  other,  as  opposed  to  being  drawn 
apart.  The  load  transfer  between  the  two  members  would  then  be 
somewhat  removed  from  the  rivet  line  and  distributed  near  the  web 
end  of  the  fleinge,  hence,  the  desirability  of  keeping  the  flange 
width  small. 

The  above  recommendations  are  In  no  way  Intended  as  reconaenda- 
tlons  for  wing  construction  as  such,  but  are  suggestions  whereby 
closer  agreement  between^analytlcal  and  experimental  results  might 
be  realized. 
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BINDING  TEST 
P   =  Dynajnoaeter  Strain  Readings.    B,  C,  S,  F,  and  G  =  Strain  Gage  Readings. 

P  B  C  D 


362.8 

689.7 

485.0 

472.7 

439.9 

417  7 

394.8 

754.9 

429.0 

242.0 

•a       36U.2 

619  6 

484.8 

472.1 

444.1 

419.5 

392.8 

754.4 

430.4 

243.8 

5       365.6 

695.2 

488.3 

471.0 

455.5 

423.7 

389.2 

733.5 

433.0 

247.6 

«  «  366.9 

697.9 

490.2 

470.0 

460.1 

426.8 

387.3 

751.6 

434.5 

249.9 

a  5  569.1 

704.9 

495.9 

472.4 

468.4 

431.0 

584.9 

750.3 

436.6 

255.2 

6  S  371.0 

707.3 

496.6 

470.7 

475.5 

434.9 

384.5 

729.6 

438.8 

255.8 

^  e^  375.2 

711-2 

499.0 

466.1 

491.1 

442.1 

377.9 

726.7 

443.2 

262.4 

t      380. 0 

720  5 

505.7 

465.8 

507.0 

449.3 

371.6 

724.9 

448.3 

269.5 

w       38iv.O 

727.5 

510  3 

463.6 

521.7 

456.8 

566.7 

722.3 

452.4 

276.2 

«  5  366.4 

691.8 

487.8 

479.8 

430.3 

413.3 

397.6 

737.7 

428.0 

239.0 

1  I  371  0 

678.8 

477.4 

476.9 

412.1 

404.2 

402.2 

739.6 

422.1 

250.8 

5  £376.0 

671  5 

472.3 

480.2 

591  5 

395.1 

408.5 

742.6 

417.2 

223.5 

5  §  386.0 

0     0 

656.9 

461.9 

484.3 

551.5 

377.4 

419.6 

747.7 

406.7 

208.9 

10 


c    562.8 

^       364.2 

684.4 

508.8 

463.0 

619.3 

562.5 

357.5 

658.3 

476.8 

452.2 

689.0 

510.3 

461.9 

625.2 

364.2 

335.9 

667.7 

479.5 

448.8 

5       565.6 

698.0 

514.0 

459.0 

636.6 

367.8 

351.6 

686.2 

485.2 

442.1 

w  ^  366.9 

704.5 

516.6 

457.1 

644.6 

370  5 

328.6 

699.2 

487.9 

436.7 

§,  5  369.1 

712.6 

520.2 

455.0 

658.3 

375.7 

^4.9 

715.7 

493.6 

429.5 

8  a  371.0 

718.8 

522.9 

454.1 

665.0 

577.3 

322.2 

728.3 

497.4 

424.8 

5  ^  375.2 

755.0 

529.8 

450.5 

683.5 

581.5 

515.5 

760.3 

504.2 

407.7 

2     380.0 

751.8 

557.2 

447.4 

705.1 

585.7 

304.1 

791.5 

510.2 

388.9 

«       384.0 

767.2 

544.3 

445.3 

725.8 

589.5 

295.2 

820.6 

515.5 

368.6 

8)       §  ^'^ 

675.5 

505.3 

466.5 

612.1 

559.2 

541.9 

638.4 

469.7 

458.2 

<S       -3  371.0 

656.7 

497.6 

471  5 

584.9 

549.9 

547.4 

600.9 

454.0 

464.6 

fl  5  S  576.0 

636.9 

490.8 

479.9 

560.2 

559.7 

554.6 

558.6 

434.0 

467.4 

1^   1  586.0 
w  w  0 

600.0 

477.0 

494.1 

516.6 

321.2 

360.2 

490.3 

394.8 

459.0 

Strain  Readings  In  Micro  Inches 


Table  I 


47 


Table  II 
BEMDIlfG  TEST 

STRAIH  QACaCS  (all  strain  readings       10     ) 


Dyn. 

lbs. 

B 

1 

1 

2 

2 

5 

3 

562.2 

0 

562.8 

12.0 

0 

698.7 

0 

485.0 

0 

472.7 

0 

56U.2 

41.0 

29.0 

691.6 

1.9 

484.8 

-.2 

472.1 

-  .6 

565.6 

69.0 

57.0 

695.2 

5.5 

488.3 

5.5 

471.0 

-1.7 

366.9 

95.0 

83.0 

697.9 

8.2 

490.2 

5.2 

470.0 

-2.7 

569.1 

lUO.O 

128.0 

704.9 

15.2 

495-9 

10.9 

472.1 

-  .6 

571  0 

178.0 

166.0 

707.5 

17.6 

496.6 

11.6 

470.7 

-2.0 

575.2 

262.0 

250.0 

711.2 

21-5 

499.0 

14.0 

466.1 

-6.6 

580.8 

575.0 

565.0 

720.5 

50.8 

505.7 

20.7 

465.8 

-6.9 

584.0 

457.0 

425.0 

727.5 

57.8 

510.5 

25.5 

465.6 

-9.1 

Dyn.       lbs. 


562.2       0 

562.8  12.0  0  439.9         0  417.7         0  394.8  0 
564.2  41.0  29.0  444.1        4.2  419.5        1.8  592.8  -2.0 
565.6  69.0  57.0  455.5  13.6  423.7        6.0  389.2  -5.6 

566.9  95.0  83.0  460.1  20.2  426.8         9.I  387.5  -7.5 

569.1  140.0  128.0  468.4  28.5  431.0  13.3  584.9  -9.9 
371.0  178.0  166.0  475.5  35.6  434.9  17.2  384.5  -10.3 

375.2  262.0  250.0  491.1  51.2  442.1  24.4  377.9  -16.9 
580.8  375.0  363.0  507.0  67.1  449.3  51.6  571.6  -23.2 
584.0  437.0  425.0  521.7  81.8  456.8  39.1  566.7  -28.1 

Dyn.  lbs.  D 


562.2 

0 

562.8 

12.0 

0 

754.9 

0 

429.9 

0 

242.0 

0 

564.2 

41.0 

29.0 

754.4 

-.5 

430.4 

1.4 

243.8 

1.8 

365.6 

69.0 

57.0 

755.5 

-1.6 

435.0 

4.0 

247.6 

5.6 

566.9 

95.0 

83.0 

751.6 

-5.5 

454.5 

5.5 

249.9 

7.9 

369.1 

140.0 

128.0 

750.5 

-4.6 

436.6 

7.6 

255.2 

11.2 

571.0 

178,0 

166.0 

729.6 

-5.5 

438.8 

9.8 

255.8 

15.8 

575.2 

262.0 

250.0 

726.7 

-8.2 

443.2 

14.2 

262.4 

20.4 

380.8 

575.0 

565.0 

724.9 

-10.0 

448.3 

19.5 

269.5 

27.5 

384.0 

457.0 

425.0 

722.3 

-12.6 

452.4 

23.4 

276.? 

54.2 
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Table  II  BIEIIDIHG  TEST  (Cont.) 


Syn.       lbs . 


562.2 

0 

562.8 

12.0 

C 

619.5 

0 

562.5 

0 

557.5 

0 

56U.2 

Ul.O 

29.0 

625.2 

5.9 

56U.2 

1.7 

555.9 

-l.k 

565.6 

69.0 

57.0 

656.6 

17.5 

567.8 

5.5 

551.6 

-5.7 

566.9 

95.0 

85.0 

skk.e 

25.5 

570.5 

8.0 

528.6 

-8.7 

569.1 

li^0.0 

128.0 

658.5 

59.  Q 

575.7 

11.2 

52H.9 

-12. if 

571.0 

178.0 

166.0 

665.0 

1+5.7 

577.5 

lk,Q 

522.2 

-15.1 

575.2 

262.0 

250.0 

685.5 

6U.0 

581.5 

19.0 

515.5 

-25.8 

580.8 

.575.0 

565.0 

705.1 

85.8 

585-7 

25.2 

50it.l 

-55.2 

58if.O 

^57.0 

425.0 

725.8 

10U.5 

589.5 

27.0 

295.2 

-lf2.1 

Dyn.       lbs. 


562.2 

0 

562.8 

12.0 

0 

658.5 

0 

U76.8 

0 

1*52.2 

0 

56U.2 

in.o 

29.0 

667.7 

9>k 

1*79.5 

2.7 

1*1*8.8 

-5.1+ 

565.6 

69.0 

57.0 

686.2 

27.9 

U85.2 

8.1* 

1*1*2.1 

-10.1 

566.9 

95.0 

85.0 

699.2 

1*0  9 

1*87.9 

11.1 

1+56.7 

-15.5 

569.1 

lUO.O 

128.0 

715.7 

57.»* 

1*95.6 

16.8 

1+29.5 

-22.7 

571- 0 

178.0 

166.0 

728.5 

70.0 

ksri.h 

20.6 

1*21*.  8 

-27.1* 

575.2 

262.0 

250.0 

760.5 

102.0 

501*.  2 

27.1* 

1*07.8 

-1*1*.  5 

580.8 

575.0 

565.0 

791-5 

155.0 

510.5 

55.5 

588.9 

-65-5 

58U.O 

U57.O 

U25.O 

820.6 

162.5 

515.5 

58.7 

568.6 

-85.6 
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TORSIOK  TIST 


P   =  Dynamoaetor  Strain  Raadlnga 

B,  C,  D,  D,  F,  and  G  =  Strain  Gage  Readings 


p 

B 

C 

D 

1 

2 

3 

1 

2 

5 

1 

2 

5 

35^.5 

691.5 

k9h,7 

U8P.7 

455.6 

416.8 

399.8 

752.9 

425-7 

257.7 

555-9 

690.1 

1+95.7 

1^80.0 

455.0 

416.7 

599.2 

755.2 

425-7 

257.7 

557.5 

690.2 

H99.5 

1+80.0 

1*55.4 

427-6 

597.5 

755.8 

417.7 

236.8 

361.7 

690.2 

512.1 

U79.5 

U56.9 

451.2 

591.5 

755.5 

405.0 

255.5 

565.6 

690.  U 

52i+.5 

i+79.0 

457.9 

474.1 

587.6 

756.8 

591.8 

234.7 

568.5 

692.5 

552.1 

W.5 

458.8 

487.5 

384.1 

757-7 

385.0 

254.2 

P 

I 

F 

G 

1 

2 

3 

1 

2 

3 

1 

2 

3 

55^.5 

680.7 

508. 1 

J+65.6 

615-5 

558.9 

357.7 

648.9 

470.9 

450.6 

555.9 

680.5 

508.7 

1^66.1 

615.9 

558.8 

557.9 

648.9 

470.3 

450.4 

557-5 

680.9 

51^.  1+ 

U61f.7 

614.5 

571.5 

557.5 

647-7 

482.5 

450.8 

561.7 

681.7 

529.0 

1*63.5 

614.5 

596.5 

555.8 

646.7 

516.5 

449.9 

565.6 

682.8 

^kk,Q 

k6^,6 

615.9 

417.2 

527.4 

645.8 

528.1 

446.6 

568.5 

685.9 

55I+.I 

k6k.k 

614.7 

427.8 

521.6 

645-7 

559.9 

445.5 

Strain  Readings  in  Microinclies  Table  III 
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CALCULATIOK  OF  EXPERIMEITTAL  SHEARHIG  STRAINS 
The  doslred  shearing  strains  ceua  be  obtained  from  the  following 
relations:  h  ^  -       z  -   — ^ — ^   vhere   i,  zt       3  =  strain  readings 


or   r  =  2   2  -  (  1+  3) 


r  =  shearing  strain 


Table  IV 

Shearing  Strains  for  Bending  Test 
V  =  Applied  load  In  pounds. 

Gage  C 

2  2  10^  -(  1+  3)  10^  r    10^  r   ,  10^ 

exp  anal 

0           0  0 

2.2        l.l^  l.Ul 

8.0         k,0  2.90 

12.7          5.5  1^.20 

18.6  8.0  6.60 
25.5  9.1  8.50 
51^.3  1^.5  12.81 
1*3.9         29.5  18.55 

53.7  21^.5  21.75 


Gage  B 

V   2  2  10^  -(   1+3)  10^  r    10^  r   ,  10^ 

^  exp  anal 

0    0         0  0  0 

29     -•'^                  1.3  -1.7  2.10 

57     6.6                   3.8  2.8  if. 20 

83    10. If        5.5  U.9  6.22 

128    21.8       llv.6  7-2  9.55 

166    23.2       15-6  7.6  12.25 

250   28.0      14.9  13.1  18.50 

363   iflA      23.9  17.5  26.90 

425   506      28.7  21.9  31-50 


0 

0 

29 

3.6 

57 

12.0 

83 

18.2 

128 

26.6 

166 

34.4 

250 

48.8 

363 

73.2 

425 

78.2 

c 


1 


. 


Gage  D 
V   2  2  10^  -(   i"*-  3)  10^  r     10^  r   ,  10^ 


•xp       anal 


0  0          0 

29  2.8  1.3 

57  8.0  k.O 

85  11  0  k.6 

128  15.2  6,6 

166  19.6  8.5 

250  28.1*  12.2 

363  38.6  17.3 

425  46.8  22. U 

Gege  G 

0  0          0 

29  5".  4  6.0 

57  16.8  17.8 

85  22.3  25.4 

128  53.6  34.7 

166  Ul.2  U2.6 

250  54.8  57.5 

363  67.0  69.7 

U25  77.4  78.7 

Gage  F 

0  0          0 

29  3.4  4.5 

57  10.6  11.6 

83  16.0  16.6 

128  22.4  26.6 

166  29.6  30- 6 

250  58.0  40.2 

363  46.4  52.6 

425  54.0  62.4 


0 

0 

1.5 

2.10 

4.0 

4.20 

5.4 

6.22 

8.6 

9.55 

11.1 

12.25 

16.2 

18.5 

21.3 

26.9 

24.4 

31.6 

0 

0 

-.6 

0 

-1.0 

0 

-3.2 

0 

-1.1 

0 

-1.4 

0 

-2.7 

0 

-2.7 

0 

-1.3 

0 

0 

0 

-1.1 

0 

-1.0 

0 

-0.$ 

0 

-4.2 

0 

-1.0 

0 

-2.2 

0 

-6.2 

0 

-8.4 

0 
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